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ABSTRACT. A new method for finding small inhomogeneities from surface scattering data is 
proposed and mathematically justified. The method allows one to find small holes and cracks 
in metallic and other obejcts from the observation of the acoustic field scattered by the objects. 



1. Introduction. 

In many applications one is interested in finding small inhomogeneities in a medium from 
the observation of the scattered field, acoustic or electromagnetic, on the surface of the 
medium. 

We have two typical examples of such problems in mind. The first one is in the area of 
material science and technology. Suppose that a piece of metal or other material is given 
and one wants to examine if it has small cavities (holes or cracks ) inside. One irradiates 
the metal by acoustic waves and observes on the surface of the metal the scattered field. 
From these data one wants to determine: 

1) are there small cavities inside the metal? 

2) if there are cavities, then where are they located and what are their sizes? 

Similar questions can be posed concerning localization not only of the cavities, but any 
small in comparison with the wavelength, inhomogeneities. Our methods allow one to answer 
such questions. 

As a second example, we mention the mammography problem. Currently x-ray mammog- 
raphy is widely used as a method of early diagnistics of breast cancer in women. However, 
it is believed that the probability for a woman to get a new cancer cell in her breast as 
a result of an x-ray mammography test is rather high (about 35 percent). Therefore it is 
quite important to introduce ultrasound mammography tests. This is being done currently. 
A new cancer cells can be considered as small inhomogeneities in the healthy breast tissue. 
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The problem is to localize them from the observation on the surface of the breast of the 
scattered acoustic field. 

The purpose of this short paper is to describe a new idea of solving the problem of 
finding inhomogeneities, small in comparison with the wavelength, from the observation of 
the scattered acoustic or electromagnetic waves on the surface of the medium. 

For simplicity we present the basic ideas in the case of acoustic wave scattering. These 
ideas are based on the earlier results on wave scattering theory by small bodies, presented 
in [l]-[5]. Our objective in solving the inverse scattering problem of finding small inhomo- 
geneities from surface scattering data are: 

1) to develop a computationally simple and stable method for a partial solution to the 
above inverse scattering problem. The exact inversion procedures (see [3] and references 
therein) are computationally difficult and unstable. In practice it is often quite important, 
and sometimes sufficient for practical purposes, to get a "partial inversion", that is, to 
answer questions of the type we asked above: given the scattering data, can one determine 
if these data correspond to some small inhomogeneities inside the body? If yes, where are 
these inhomogeneities located? What are their intensities? We define the notion of intensity 
v m of an inhomogeneity below formula (1). 

In section 2 the basic idea of our approach is described. In section 3 its short justification 
is presented. Some theoretical and numerical results based on a version of the proposed 
approach one can find in [5] . 

2. Basic equations. 

Let the governing equation be 

[V 2 + k 2 + k 2 v(x)]u = -5(x-y) in M 3 , (1) 

where u satisfies the radiation condition, k = const > 0, and v(x) is the inhomogeneity in 
the velocity profile. 

Assume that sup^^a \v(x)\ < c , supp v = U^ =1 B m (z m , p m ) C R 3 _ = {x : X( 3 ) < 0}, 
where X( 3 ) denotes the third component of vector x in Cartesian coordinates, B m (z m , p m ) 
is a ball, centered at z m with radius p m , kp m <C 1. 

Denote 

:= / v{x)dx. 

Inverse Problem (IP). Given u(x, y, k) for all x, y G P, P = {x : X( 3 ) = 0} and a fixed k > 0, 
find {z m ,v m }, 1 < m < M. 

In this paper we propose a numerical method for solving the (IP). 
To describe this method let us introduce the following notations: 

P := {x : x (3) = 0}, (2) 

{xj,yj} := £j, l<j< J, xj,yj G P, 

are the points at which the data u(xj,yj, k) are collected, (3) 
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k > is fixed, 



(4) 



g(x,y,k) :-- 



exp(ik\x — y\) 
Air\x — y\ 



(5) 



Gj(z) :=G(^,z) := g(x j ,z,k)g(y j ,z,k), 



f.i 



u(x j ,y j ,k)-g(x j ,y j ,k) 
k 2 



(6) 



(7) 



$(zi,. ..,z m , Wi, 



M 



m=l 



(8) 



The proposed method for solving the (IP) consists in finding the global minimizer of 
function (8). This minimizer (z\, . . . , z m , t>i, . . . ,v m ) gives the estimates of the positions 
z m of the small inhomogeneities and their intensities v m . This is explained in more detail 
below formula (14). Numerical realization of the proposed method, including a numerical 
procedure for estimating the number M of small inhomogeneities from the surface scattering 
data is described in [8] . 

Our approach with a suitable modification is valid in the situation when the Born ap- 
proximation fails, for example, in the case of scattering by delta-type inhomogeneities [9]. 

In this case the basic condition Mk 2 cop 2 « 1 (*) which guarantees the applicability of 
the Born approximation is violated. Here p := maxi< m <M p m and Co was defined below 
formula (1). We assume throughout that M is not very large, between 1 and 15. 

In the scattering by a delta- type inhomogeneity the assumption is cop 3 = const := V as 
p — > 0, so that for any fixed k > one has k 2 cop 2 = k 2 V p~ x — > oo as p — » 0, and clearly 
condition (*) is violated. 

In our notations this delta-type inhomogeneity is of the form k 2 v(x) = k 2 J2m=i ^mH x ~ 

The scattering theory by the delta- type potentials (see [9]) requires some facts from the 
theory of selfadjoint extensions of symmetric operators in Hilbert spaces and in this short 
paper we will not go into detail. 

3. Justification of the proposed method. 



We start with an exact integral equation equivalent to equation (1) with the radiation 
condition: 

M 

(9) 

m= 



M ,. 

u(x,y,k) = g(x,y,k) + k 2 ^2 / g(x,z,k)v(z)u(z,y,k)dz. 

m=l 



For small inhomogeneities the integral on the right-hand side of (9) can be approximately 
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written as 



g(x,z,k)v(z)u(z,y,k)dz := k 2 / g(x, z,k)v(z)g(z,y,k)dz + e 2 (10) 

J Bm 

= k 2 G(x,y,z m ) / vdz + 



e 2 



= k z G(£,z m )v m + e 2 , l<m<M, 

where e 2 is denned by the first equation in formula (10), it is the error due to replacing u 
under the sign of integral in (9) by g, ~z m is a point close to z m . 

One has \u— g\ = 0( Mk d %" p ) ifx,y G P, so the error term e 2 in (10) equals to 0( M k d s oP ) 
if x, y G P. 

Therefore the function u(z, y, k) under the sign of the integral in (9) can be replaced by 
g(x, y, k) with a small error provided that 

k 2 p 3 

c M— p<l, x,yeP, (11) 

where p = maxi< m <M p m , c := max xeR 3 |f M is the number of inhomogeneities, d 
is the minimal distance from B m , m = 1, 2, M to the surface P, d ~ 1 means that the 
length is measured in the units of length d. 

If a sufficient condition for the validity of the Born approximation holds, that is, 

Mk 2 c p 2 := 6 « 1, (12) 

then 

JM-2L.4 2 6 x2 2 

0(^^-)=0(^-)«l ifd~l. 

Note that u in (9) has dimension L _1 , where L is the length. 

If the Born approximation is not valid, for example, if cop 3 = V ^ as p — > 0, which is 
the case of scattering by delta-type inhomogeneities, then the error term e 2 in formula (10) 

2 4 2 2 4 2 

can still be negligible: in this case e 2 = 0( M % 3 V ), so e 2 « 1 if M % v « 1. 

If one understands a sufficient condition for the validity of the Born approximation as 
the condition which guarantees the smallness of e 2 for all x,y G M 3 then condition (12) is 
such a condition. However, if one understands a sufficient condition for the validity of the 
Born approximation as the condition which guarantees the smallness of e 2 for x, y running 
only through the region where the scattered field is measured, in our case when x,y G P, 
then a much weaker condition (11) will suffice. 

In the limit p — » and cop 3 = V ^ formula (10) takes the form (13) (see below). This 
can be derived from [9, p. 113]. Formula (1.1.33) in [9] shows that the resolvent kernel of 
the Schrodinger operator with the delta-type potential supported on a finite set of points 
(in our case on the set of points Z\, ....,zm) has the form 

M 

u(x, y, k) = g(x, y, k) + k 2 ^ c mm >g{x, z m )g(y, z m >), (13) 

m=l 
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where c mm / are some constants. These constants are determined by a selfadjoint real- 
ization of the corresponding Schrodinger operator with delta-type potential. There is an 
M 2 — parametric family of such realizations (see [9] for more details). 

Although in general the matrix c mm < is not diagonal, under a practically reasonable 
assumption (11) one can neglect the off-diagonal terms of the matrix c mm i and then formula 
(13) reduces practically to the form (10) with the term e 2 neglected. 

We have assumed in (10) that the point z^ exists such that f B g(x, z, k)v(z)g(z, y, k)dz = 
G(x, y, z^)u^. This is an equation of the type of mean-value theorem. However, such a theo- 
rem does not hold, in general, for complex-valued functions. Therefore, if one wishes to have 
a rigorous derivation, one has to add to the error term e 2 in (10) the error which comes from 
replacing of the integral f B g(x, z, k)v(z)g(z, y, k)dz in (10) by the term G(x, y, z^)v^. The 
error of such an approximation can be easily estimated. We do not give such an estimate, 
because the basic conclusion that the error term is negligible compared with the main term 
k 2 G(x,y,z^)u^ remains valid under our basic assumption kp << 1. From (10) and (7) it 
follows that 

M 

fj ~ ^2 Gj{z m )v m , Gj(z m ) := G(£j,z m ,k) (14) 

771 = 1 

Therefore, parameters z m and v m can be estimated by the least-squares method if one finds 
the global minimum of the function (8): 

$(zi, . . . ,z M , vi, . . . ,v M ) = min. (15) 

Indeed, if one neglects the error of the approximation (10), then the function (8) is a smooth 
function of several variables, namely, of z\, z%, ■■■■Z m , Vi, V2, —v m , and the global minimum 

of this function is zero and is attained at the actual intensities Vi,V2, ,v m and at the 

values Zi = z~i,i = 1, 2, ....to,. 

This follows from the simple argument: if the error of approximation is neglected, then the 
approximate equality in (14) becomes an exact one. Therefore fj — X)m=i Gj(z m )v m = 0, so 
that function (8) equals to zero. Since this function is non-negative by definition, it follows 
that the values ~z m and v m are global minimizers of the function (8). Therefore we take the 
global minimizers of function (8) as approximate values of the positions and intensities of 
the small inhomogeneities. 

In general we do not know that the global minimizer is unique. For the case of one small 
inhomogeneity (to = 1) uniqueness of the global minimizer is proved in [5] for all sufficiently 
small p m for a problem with a different functional. The problem considered in [5] is the 
(IP) with M = 1, and the functional minimized in [5] is specific for one inhomogeneity. 

The scattering theory for small scatterers originated in the classical works of Lord Rayleigh. 
It was developed in [1] and [2], where analytical formulas for the scattering matrix were 
derived for the acoustic and electromagnetic scattering problems. In [1] and [3] inverse scat- 
tering problems for small bodies are considered. Numerically an important ingredient of 
our approach is the solution of the global minimization problem (14). The theory of global 
minimization is developed extensively and the literature of this subject is quite large. We 
mention two recent papers [6] and [7] where the reader can find many references. 
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